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In this paper we consider three dimensional BTZ black holes with three models of nonlinear 
electrodynamics as source. Calculating heat capacity, we study the stability and phase transitions 
of these black holes. We show that Maxwell, logarithmic and exponential theories yield only type 
one phase transition which is related to the root(s) of heat capacity. Whereas for correction form 
of nonlinear electrodynamics, heat capacity contains two roots and one divergence point. Next, we 
use geometrical approach for studying classical thermodynamical behavior of the system. We show 
that Weinhold and Ruppeiner metrics fail to provide fruitful results and the consequences of the 
Quevedo approach are not completely matched to the heat capacity results. Then, we employ a new 
metric for solving this problem. We show that this approach is successful and all divergencies of its 
Ricci scalar and phase transition points coincide. We also show that there is no phase transition for 
uncharged BTZ black holes. 


I. INTRODUCTION 

One of the interesting subjects for recent gravitational studies is the investigation of three dimensional black holes 
Considering three dimensional solutions helps us to find a profound insight in the black hole physics, quantum 
view of gravity and also its relations to string theory [iMl- Moreover, three dimensional spacetimes play an essential 
role to improve our understanding of gravitational interaction in low dimensional manifolds [T5j |. Due to these facts, 
some of physicists have an interest in the (2 + l)-dimensional manifolds and their attractive properties [TgI - I?^ . 

The Maxwell theory is in agreement with experimental results, but it fails regarding some important issues such as 
self energy of point-like charges which motivates one to regard nonlinear electrodynamics (NED). There are some evi¬ 
dences that motivate one to consider NED theories: solving the problem of point-like charge self energy, understanding 
the nature of different complex systems, obtaining more information and insight regarding to quantum gravity, com¬ 
patible with AdS/CFT correspondence and string theory frames, description of pair creation for Hawking radiation 
and the behavior of the compact astrophysical objects such as neutron stars and pulsars (28l - [30l| . Therefore, many 
authors investigated the black hole solutions with nonlinear sources [sTI - ld^ . 

On the other hand, thermodynamical structure of the black holes, has been of great interest. It is due to the fact 
that, according to AdS/CFT correspondence, black hole thermodynamics provides a machinery to map a solution in 
AdS spacetime to a conformal field on the boundary of this spacetime |13l - [l5l| . Also, it was recently pointed out that 
considering cosmological constant as a thermodynamical variable leads to the behavior similar to the Van Der Waals 
liquid/gas system |50l - l57j| . In addition, phase transition of the black holes plays an important role in exploring the 
critical behavior of the system near critical points. There are several approaches that one can employ to study the 
phase transition. One of these approaches is studying the behavior of the heat capacity It is argued that roots and 
divergence points of the heat capacity are representing two types of phase transition (58l - l64| . In addition, studying 
the heat capacity and its behavior, enable one to study the thermal stability of the black holes . 

Another approach for studying phase transition of black holes is through thermodynamical geometry. The concept 
is to construct a spacetime by employing the thermodynamical properties of the system. Then, by studying the 
divergence points of thermodynamical Ricci scalar of the metric, one can investigate phase transition points. In 
other words, it is expected that divergencies of thermodynamical Ricci scalar (TRS) coincide with phase transition 
points of the black holes. Firstly, Weinhold introduced differential geometric concepts into ordinary thermodynamics 
[ 71 L l72j | . He considered a kind of metric defined as the second derivatives of internal energy with respect to entropy 
and other extensive quantities for a thermodynamical system. Later Ruppeiner [zlizi introduced another metric 
and defined the minus second derivatives of entropy with respect to the internal energy and other extensive quantities. 
The Ruppeiner metric is conformal to the Weinhold metric with the inverse temperature as the conformal factor. It 
is notable that, both metrics have been applied to study the thermodynamical geometry of ordinary systems (75l - [80l| . 
In particular, it was found that the Ruppeiner geometry carries information of phase structure of thermodynamical 
system. For the systems with no statistical mechanical interactions (for example, ideal gas), the scalar curvature 
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is zero and the Ruppeiner metric is flat. Because of the success of their applications to ordinary thermodynamical 
systems, they have also been used to study black hole phase structures and lots of results have been obtained for 
different types of black holes [sTI - l^ . 

It is notable that these two approaches fail in order to describe phase transition of several black holes 87|. In order 
to overcome this problem, Quevedo proposed new types of metrics for studying geometrical structure of the black hole 


thermodynamics 
the black holes [d 
But this a 
capacity [8 


This method was employed to study the geometrical structure of the phase transition of 
494 1 and proved to be a strong machinery for describing phase transition of the black holes. 


roach was not completely coincided with the results of classical thermodynamics arisen from the heat 
In Ref. [s^ , a new metric was proposed in which the denominator of its Ricci scalar is only constructed 
of numerator and denominator of the heat capacity. Several phase transition of the black holes have been studied in 
context of the HPEM (Hendi-Panahiyan-Eslam Panah-Momennia) metric [s^ . 

In this paper we study thermal stability and phase transition of the BTZ black holes in presence of several NED 
models in context of heat capacity. Then, we employ Weinhold, Ruppeiner and Quevedo methods for studying 
geometrothermodynamics of these black holes. We will see that Weinhold and Ruppeiner metrics fail to provide 
fruitful results and the consequences of the Quevedo approach are not completely matched to the heat capacity 
results. Then, we employ the HPEM metric and study the phase transition of these black holes in context of 
geometrothermodynamics. We end the paper with some closing remarks. 


II. NONLINEARLY CHARGED BTZ BLACK HOLE SOLUTIONS 

The (2 + l)-dimensional action of Einstein gravity with NED field in the presence of cosmological constant is given 
by 


/ = - y d^x^g [i? - 2A + L{T )], 


( 1 ) 


where R is the Ricci scalar, the cosmological constant is A = —l~^ in which I is a scale factor. Also, 
Lagrangian of NED, in which we consider three models. First model was proposed by Hendi (HNED) 
one is Soleng theory (SNED) [9^ and third one is correction form of NED (CNED) 


L(F) is the 
[ 93 , second 


L{F) = 


/32 (exp(-^) - 1 ), 

-8/3" In (1 + ^) , 


-F + aF^ + 0{a^), 


HNED 

SNED > 
CNED 


( 2 ) 


where /3 and a are called the nonlinearity parameters, the Maxwell invariant F = in which = d^Ay — dyA^ 

is the electromagnetic field tensor and A^ is the gauge potential. We should note that for [3 —>■ 00 (HNED and SNED 
branches) and a —> 0 (CNED branch) the Maxwell Lagrangian can be recovered. 

The nonlinearly charged static black hole solutions can be introduced with the following line element 

ds^ = -f{r)df + + r'^de'^, (3) 

fir) 

where the metric function /(r) was obtained in Refs. [zaiii 


= ]2 - ™ + ' 
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( 4 ) 


where P = yl + m and q are integrations constant which are related to mass parameter and the electric charge 
of the black hole, respectively. In addition, Lw = LambertW 7 = 7 ( 0 ) ~ 0.57722 and the special function 

Ei{l,x) = ^dz. 
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The entropy and the electric charge of the obtained NED black hole solutions can be calculated with the following 
forms [^. 1^ 


S = 
Q = 


7rr+ 

TTq 

T’ 


(5) 

( 6 ) 


where r+ is the event horizon of black hole. On the other hand, the quasi-local mass, which is related to geometrical 
mass, can be obtained as Hil 


M = —m. 


(7) 


Regarding Eqs. dSD and ([6]) with ([7]) and obtaining m by using /(r = r+) = 0, for these three cases of BTZ black 
holes we can write 


4/3SQ(1-2L'^) 2S=(/^/3^-2) 


M = - X < 

Stt 


45"^ 


16/3^5*2 In + 3 


p --+4Q2 

' -4Q2 
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( 8 ) 


where E' = \J^ + and = LambertW . Having conserved and thermodynamic quantities at hand, it 

was shown that the first law of thermodynamics may be satisfied [9^. The main goal of this paper is investigating 
thermal stability and phase transition for these black holes. 


A. Heat capacity 


In order to investigate thermal stability and phase transition, one can usually adopt two different approaches to the 
matter at hand. In one method, the electric charge is considered as a fixed parameter and heat capacity of the black 
hole will be calculated. The positivity of the heat capacity is sufficient to ensure the local thermal stability of the 
solutions and its divergencies are corresponding to the phase transition points. This approach is known as canonical 
ensemble. Another approach for studying thermal stability of the black holes is grand canonical ensemble. In this 
approach, thermal stability is investigated by calculating the determinant of Hessian matrix of M(S,Q) with respect to 
its extensive variables. The positivity of this determinant also represents the local stability of the solutions. Although 
these two approaches are different fundamentally, one expects that the results being the same for both ensembles; i.e. 
ensemble independent. Here we use the first method for studying thermal stability. For this purpose, the system is 
considered to be in fixed charge and the heat capacity has the following form 



-1 

Q 


(9) 


It is notable that, when we study heat capacity for investigating the phase transition, we encounter with two different 
phenomena. In one, the changes in the signature of the heat capacity is representing a phase transition of the system. 
In other words, if the heat capacity is negative, then the system is in thermally unstable phase, whereas for the case 
of the positive Cq, the system is thermally stable. The roots of the heat capacity in this case are representing phase 
transition points which means one should solve the following equation 



dM 
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( 10 ) 
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where hereafter we call this type of the phase transition as type one. It is a matter of calculation to show that 

2nSPy/f^{l + L{y) 


S^/3='(l-S^VF)-6Q*i^ 25=^ l„(r:l^)+6/3=+Q=' 


SNED 


{2Q*a-S^Q'^)l^+S* 


CNED 


The other case of phase transition is the divergency of the heat capacity. In other words, the singular points of 
the heat capacity are representing places in which system goes under phase transition. This assumption leads to the 
fact that the roots of the denominator of the heat capacity are representing phase transitions. Therefore, we have the 
following relation for this type of phase transition 

where we call this type of the phase transition as type two. Due to economical reasons, we did not write the explicit 


relations of 




for different BTZ solutions. 


B. Geometrical Thermodynamics 

In order to have an effective geometrical approach for studying phase transition of a system, one can build a suitable 
thermodynamical metric and investigate its Ricci scalar. Thermodynamical metrics were introduced based on the 
Hessian matrix of the mass (internal energy) with respect to the extensive variables. Therefore, although the electric 
charge is a fixed parameter for calculating the heat capacity in canonical ensemble, it may be an extensive variable for 
constructing thermodynamical metrics. In this method we expect that TRS diverges in both types of the mentioned 
phase transition points. In other words, the denominator of TRS must be constructed in a way that contains roots of 
the denominator and numerator of the heat capacity. In what follows, we will study the denominator of TRS of the 
several geometrical approaches and follow the recently proposed thermodynamical metric which its denominator only 
contains numerator and denominator of the heat capacity, and therefore, divergencies of TRS coincide with roots and 
divergences of the heat capacity. 

In order to find the roots and divergence points of heat capacity, we should solve its numerator and denominator, 
separately. Solving the mentioned equations with respect to entropy, leads to 


^0 ^ = 
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—2 exp 
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It is evident from obtained equation for HNED and SNED that there is only one real positive entropy in which 
heat capacity vanishes. Interestingly, in case of CNED theory we find two roots for heat capacity. It is evident that 
the roots are increasing functions of the electric charge in these theories. As for nonlinearity parameter, in case of the 
HNED and SNED theories, the root is an increasing function of /3. Whereas for CNED theory, the smaller root is an 
increasing function of a whereas the larger root is a decreasing function of it. 

Now we are in position to study the existence of the type two phase transition point which is related to divergency of 
the heat capacity. Considering HNED branch of Eq. (HI, one finds Soo is not real for all values of I and /3. Therefore, 
there is no physical divergence point for the heat capacity of HNED model. Next, for the case of SNED model the 
same behavior is seen. In other words, since 0 < exp(^) < I for 0 < a: < 00 , we can not obtain real Sea for all values 
of I and p. Next, we should investigate CNED model. Regarding Eq. da, we find that there is a divergence point 
for heat capacity (we should note that in this paper we consider positive a). In other words, this theory of nonlinear 
electromagnetic field enjoys the phase transition of type two. The divergence point is an increasing function of the 
nonlinearity and electric charge parameters. It is worthwhile to mention that for the case of vanishing nonlinearity 
parameter in this theory. Sea goes to zero. In other words, there is no divergence point for heat capacity of Maxwell 
theory. It is also clear that in case of chargeless BTZ black holes, there is no root and divergence point for heat 
capacity. In other words, in case of chargeless BTZ black holes, there is not any kind of phase transition. 


1. Weinhold and Ruppeiner metrics 
The Weinhold metric was given in [zHil 

dS^ = gYbdX^^dX^ (17) 

where = d'^M (A'^) jdX^'dX^ and also A“ = A“ (S', W), where W denotes other extensive variables of the 
system. In case of Weinhold approach, one is considering the mass of the system as potential, other parameters 
such as entropy and electric charge as extensive parameters and related quantities such as temperature and electric 
potential as intensive parameters. 

The Ruppeiner metric was defined as 


dSji = g^,dX<^dX\ (18) 

where = —d^S (A'^) /9A“9A*' and A“ = A“ (M, A*). In this case the thermodynamical potential is entropy. It 
is worthwhile to mention that according to the proposal of the Quevedo, these two approaches are related to each 
other by a Legendre transformation [^ . 

Taking into account thermodynamical metrics of Weinhold and Ruppeiner, one can obtain their Ricci scalars. Since 
we would like to investigate divergence points of TRS, 72., we focus on its denominator (7772). One finds 

D^^W^n ^ M^MssMqq-MIq)\ (19) 

DTZ^^P = TM^MssMqq-MIq)\ ( 20 ) 

where Mk = ^ and Mkj = 

2. The Quevedo metrics 

The Quevedo metrics have two kinds with the following forms [H, 

dSj, = g2,dX<^dX\ 


where is 


( 21 ) 
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FIG. 1: Weinhold Ricci scalar (solid line), heat capacity (dotted line) and temperature (dashed line) versus S for I = 1. 
HNED model: (up-left panel) and SNED model: (up-right panel): q = 0.3, /3 = 1. 

CNED model: (down-left panel) and CNED model: (down-right panel): q — 1, a = 0.007 (different scales). 


with 


T = 


SMs + QMq, 

SMs, 


Case I 
Case II 


Taking into account Quevedo metrics, one can find that their related (denomerator of) Ricci scalars can be written as 

DTZ^-^ = MlsM^QiSMs + QMgf, (23) 

Dt^Q-ii = S^MlMlsMhn. (24) 


3. HPEM metric 


In order to avoid any extra divergencies in TRS which may not coincide with phase transitions of the type one and 
two, and also ensure that all the divergencies of the TRS coincide with phase transition points of the both types, 
HPEM metric was introduced [83 
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Ms ( -Mss 0 
0 M. 


Mqq 


QQ 


(25) 































7 




FIG. 2: Ruppeiner Ricci scalar (solid line), heat capacity (dotted line) and temperature (dashed line) versus S for I = 1. 
HNED model: (up-left panel) and SNED model: (up-right panel): q = 0.3, /3 = 1. 

CNED model: (down-left panel) and CNED model: (down-right panel): q — 1, a = 0.007 (different scales). 


In this case we have considered the total mass as thermodynamical potential, entropy and electric charge as extensive 
parameters. Calculations show that denominator of TRS leads to 

D-rHpem = S^MlMls. ( 26 ) 

III. THE RESULTS OF VARIOUS APPROACHES 

Here, we investigate phase transitions of black holes using geometrothermodynamics. For this purpose, we used 
thermodynamical metrics introduced in previous section for the black holes solutions obtained in the section HIl 

For Weinhold metric, none of divergencies of the Ricci scalar coincide with roots of the heat capacity in every 
theories of NED models that we have considered in this paper. On the other hand, one of the divergencies of TRS 
and divergence point of the heat capacity in CNED theory, coincide with each other. It is notable that, in cases of 
the HNED and SNED theories, there is one divergence point for TRS (up panels of Figs. [IJ whereas for CNED, there 
are two divergencies (down panels of Figs. [J). 

In case of Ruppeiner metric, for HNED and SNED models, there is a root for heat capacity in which Ricci scalar 
of the Ruppeiner metric has a divergency. But there is also another divergence point for Ricci scalar which does not 
coincide with any phase transition point (up panels of Eigs. [5]). Therefore, there is an extra divergence point. In 













































FIG. 3: Quevedo Ricci scalar case I (solid line), heat capacity (dotted line) and temperature (dashed line) versus S for I = 1. 
HNED model: (up-left panel) and SNED model: (up-right panel): q = 0.3, /3 = 1. 

CNED model: (down-left panel) and CNED model: (down-right panel): q — 1, a = 0.007 (different scales). 


case of the CNED model, two roots and one divergence point are observed for heat capacity in which Ricci scalar has 
related divergencies. In addition to these divergence points, one extra divergence point is also observed which is not 
related to any phase transition point (down panels of Figs. [2|). 

As for Quevedo metrics, for case I, similar behavior as Weinhold is observed for all three theories of NED (Fig. 

On the other hand, in case of the other metric of Quevedo, two divergence points for Ricci scalar were observed 
for HNED and SNED theories. One of these divergence points coincides with root of the heat capacity for these two 
nonlinear theories whereas the other one does not (up panels of Figs. 0]). In case of CNED theory, all the divergence 
points of TRS coincide with phase transition points except one (down panels of Figs. U]). In other words, Quevedo’s 
metric predicts an extra divergence point, corresponding to the equation Mqq = 0, which is not predicted in classical 
black hole thermodynamics. 

It is evident that in case of HPEM, all types of phase transition points of heat capacity coincide with divergencies of 
TRS of HPEM method (Figs. [S]). In other words, independent of the nonlinear theory under consideration, the HPEM 
method provide a machinery in which no extra divergency for TRS is observed and divergence points of TRS and 
phase transition points coincide. Another interesting and important property of the HPEM method is the behavior 
of TRS near divergence point for different types of phase transition. As one can see, the signature and behavior of 
TRS near divergence point for phase transition type one and two are different. Therefore, independent of studying 
the heat capacity, one can distinguish these types of phase transition from one another only by studying the behavior 
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FIG. 4: Quevedo Ricci scalar case II (solid line), heat capacity (dotted line) and temperature (dashed line) versus S for Z = 1. 
HNED model: (up-left panel) and SNED model: (up-right panel): q = 0.3, /3 = 1. 

CNED model: (down-left panel) and CNED model: (down-right panel): q = I, a = 0.007 (different scales). 

of TRS. 


IV. CLOSING REMARKS 

In this paper, we have considered BTZ black holes, in presence of three models of NED. We studied stability and 
phase transitions related to the heat capacity of the mentioned black holes. Next, we employed geometrical approach 
to study the thermodynamical behavior of the system. In other words, we have studied phase transitions of the 
system through Weinhold, Ruppeiner and Quevedo methods. Also, we used the recently proposed approach to study 
geometrical thermodynamics. 

We found that the Weinhold and Ruppeiner metrics for studying these BTZ solutions fail to provide a suitable 
result. In addition, the divergence points of the Quevedo TRS were not completely matched with the phase transition 
points of the heat capacity results. In other words, in these approaches, the existence of extra divergencies were 
observed which were not related to any phase transition point in the classical thermodynamics. In some of these 
approaches no divergency of TRS coincided with phase transition points. In order to obtain a consistent results with 
the classical thermodynamic consequences (the heat capacity), we employed a new thermodynamical metric. In this 
approach, all the divergencies of TRS coincided with phase transition points. In other words, roots and divergence 
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FIG. 5: HPEM Ricci scalar (solid line), heat capacity (dotted line) and temperature (dashed line) versus S for I = 1. 
HNED model: (up-left panel) and SNED model: (up-right panel): q = 0.3, /3 = 1. 

CNED model: (down-left panel) and CNED model: (down-right panel): q — 1, a = 0.007 (different scales). 


points of the heat capacity of the BTZ black holes in presence of each nonlinear models matched with divergencies of 
TRS of this metric. 

Also, we found that in case of HNED and SNED theories, there is no divergency for heat capacity. It means 
that, like Maxwell theory, these two theories have no second type phase transition. These two nonlinear theories 
of electrodynamics, preserved the characteristic behavior of the Maxwell theory in case of heat capacity. On other 
hand, for CNED model, two roots and one divergence point was found for the heat capacity. In other words, due to 
contribution of the nonlinear electromagnetic field, heat capacity enjoys the existence of one more phase transition 
point of type one and a phase transition of the type two. In essence this theory is a generalization of the Maxwell 
theory. But this generalization added another property to heat capacity that was not observed for the Maxwell theory. 

Finally it is worthwhile to mention a comment related to Legendre invariancy. It was shown that 98| the Leg¬ 


endre invariance alone is not sufficient to guarantee a unique description of thermodynamical metrics in terms of 
their curvatures. In addition to Legendre invariancy, one needs to demand curvature invariancy under a change of 
representation. Therefore, it will be worthwhile to investigate both Legendre and curvature invariancies. In addition, 
it will be interesting to think about the fundamental relation between the following two issues: (I) Agreement of ther¬ 
modynamical curvature results with usual thermodynamical approaches (such as the heat capacity); (II) Curvature 
invariancy in addition to the Legendre invariancy. It is also worthwhile to probe the fundamentality of cases (I) and 
(II) to find that considering which one may leads to satisfy another one. Although first issue has been investigated 
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for special cases [9^, the second one has been remained unanalyzed yet. We may address them in an independent 
work in the future. 
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